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Abstract. We calculate energies and hyperfine structure constants A and B for 
low-lying levels of 173 Yb. The latter is treated as a two-electron atom. To account 
for valence-valence correlations the configuration interaction method is used. 
Valence-core correlations are accounted for within the many-body perturbation 
theory which is used to construct effective two-electron operators in the valence 
space. 



1. Introduction 

Two years ago a method for calculation of the low-lying energy levels of many- 
electron atoms was proposed [Q. Calculations for Tl [|l], Ca, Sr. Ba and Yb Q 
and Ba || confirmed its effectiveness. Later this method was generalized for other 
observables, such as hyperfine structure (HFS) constants H and polarizabilities ||. 
Within this method valence-valence correlations are treated non-perturbatively, while 
core- valence and core-core correlations are treated within the many-body perturbation 
theory (MBPT). 

Here we report calculations for the HFS constants of the low-lying levels of 
ytterbium. Our interest to this atom is caused in part by an extremely large parity 
nonconserving (PNC) amplitude for the 1 5o (6s 2 ) — > 3 Di (5d6s) transition. It was 
first suggested by DeMille |(| and later confirmed in @, g|, that this amplitude is 
approximately 100 times larger than PNC amplitude for the 6s — > 7s transition in Cs 
and 10 times larger than that for the 6p — > 7p transition in Tl. This makes Yb a very 
promising candidate for the future PNC experiment. Therefore, an accurate atomic 
calculations for Yb are needed. Moreover, it is important to have a reliable estimate 
of the accuracy of such calculations. PNC amplitude is very sensitive to the behavior 
of the wave function in the vicinity of the nucleus. Of course, the same applies to HFS 
constants. That makes HFS calculations very important for the future analysis of the 
theoretical accuracy for the PNC amplitude. 

We also have a more general interest in ytterbium as an atom with the relatively 
shallow 4/ core shell. For such an atom it is not obvious at all, that core-valence 
correlations can be treated by means of the MBPT. In the earlier paper j!| we have 
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shown that the method works for the low-lying energy levels and here we extend our 
analysis to the HFS constants. We also show here, that with some caution it is possible 
to calculate HFS constants of the levels, which lie above the core excitation threshold. 

In the next section of the paper we define effective operators for the valence 
electrons and briefly describe how MBPT can be used to calculate these operators 
(for more details see @, |). In section | we give some details of the calculations 
followed by the discussion of the results in section 0. 



2. Effective operators for valence electrons 

At present there are several methods of calculations for many-electron atoms. 
For atoms with one electron above a closed core a conventional MBPT is quite 
effective (see, for example, ||). For atoms with several valence electrons the multi 
configurational Hartree-Fock (MCHF) method is often used |1Q| . Lately the coupled- 
claster method became very popular jOJ |l2|, All these methods were used to 

calculate HFS constants of different atoms (see, for example, Jl4|, |l5|, |l6|, [l?], 

The most complicated problem in precise atomic calculations is associated 
with the necessity to treat valence-valence correlations, core-valence, and core-core 
correlations. The former are usually too strong to be treated perturbatively, while 
other types of correlations can not be effectively treated within non-perturbative 
techniques, such as MCHF or CI method. 

Therefore, it is natural to try to combine MBPT with one of the non-perturbative 
methods. In |l[] it was suggested to use MBPT in order to construct an effective 
Hamiltonian for valence electrons. After that, Schrodinger equation for valence 
electrons is solved within CI approach. That allows to find the low-lying energy levels. 
In order to calculate other atomic observables, one need to construct corresponding 
effective operators for valence electrons Q . The latter effectively account for the core- 
valence and core-core correlation, which are of particular importance for such singular 
operators as HFS ones. 

Suppose, that we are interested in atomic levels with energies Ei — Eo < s, where 
Eq is the ground state energy. In the first approximation we can assume that inner 
electrons with the Hartree-Fock energies \e n \ ^> e form the core, which is described 
by the wave function 

tfcoro = (^ c !)" 1/2 det(0i,0 2 ,...0 JVc ), (1) 

hDF(t>n = En^n, (2) 

where /idf is the Dirac-Fock operator, and N c is the number of core electrons. Note, 
that this operator can include the field of some valence electrons as well. For example, 
it is natural to consider Yb as a two-electron atom with the core [Is 2 . . . 4/ 14 ], 
and operator /idf can be calculated for configuration Is 2 . . . 4/ 14 6s 2 (so called, V - 
approximation, N is the number of electrons in the atom). 

Let us say that many-electron wave function belongs to the valence subspace if 
core electrons are in the state (|l|). We will denote projector operator on this subspace 
by P. Then, operator Q = 1 — P will project on the subspace for which at least one 
of the core electrons is excited to one of the states 4> n with n > N c . 

One can show [Q, that Schrodinger equation 

i?* = E^ (3) 
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is equivalent to the following equation in the P-subspace for the function <J> = P\I>: 

(PiJP + = (4) 

E(P) = PV'R Q (E)V'P, (5) 

where V' is the operator of the residual Coulomb interaction and Rq(E) is the Green's 
function in Q-subspace: 

V = H - Ho, (6) 
Operator H is defined in terms of one-electron operator (||) : 

AT 

Ho = Y, h Mn)-W. (8) 

i=l 

The constant W in the right hand side of this equation is introduced to compensate 
the double counting of the two-electron interaction in the sum. It can be fixed, for 
example, by the requirement: 

(*core|-ffo|*core) = (*corc|P|*corc), (9) 

which implies, that 

N c 

W = Wo = Y^Ei - (^corclPl^corc). (10) 
i=l 

It is also possible to use W as a free parameter to match the energy spectrum. This 
subject is discussed in more details elsewhere 
Equations (|) — (0) yield: 

* = {P + R Q {E)V'P)<S>, (11) 

The orthonormality condition (^il^fc) = 8ik is approximately equivalent to the 
following condition for $: 

(*i|i-a B E(J)|* fc ) «f iit , (12) 

where E w (Ei + E k )/2. Note, that only the last of equations (Q) — (ll^) is 
approximate. 

The operator in the left hand side of equation (H) plays the role of an 
effective Hamiltonian P c ff- Equations (|J) — (fi]) allow to construct it with the 
help of conventional diagrammatic technique [Q, which is based on the following 
representation of the Green's function: 

R Q (E)=R Q (E)+R° Q (E)V'R Q (E), (13) 

R ^ E) - Q E~h^Q Q - (14) 

In particular, if we restrict ourselves to the second order MBPT, we have to replace 
Rq(E) in the right hand side of §) by R° Q {E). 

Suppose now, that we know solutions of Equation (|J) and we want to use them 
to calculate an observable a, which is associated with one-electron operator A: 

a=($\A\V}. (15) 

Let us define effective operator A e g so, that 

a=<*|Aff|*>- (16) 
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The following approximation for A c g was suggested in [Q 

A cS « PAP + PV'R° Q (E)A rpa P + PA rp& R° Q (E)V'P, (17) 

where A rpa corresponds to the well-known random-phase approximation (RPA). 

Note, that only the first order MBPT corrections are completely included in (|l7|). 
Some second order corrections, the, so-called, structural radiation, as well as many 
higher-order corrections, are omitted. Nevertheless, this approximation accounts for 
several most important MBPT corrections to all orders. Diagrammatic representation 
of Equation (|l^) can be found in [[l| . 

3. Hyperfine structure calculations 

We use basis representation for diagrams for effective Hamiltonian and for effective 
HFS operators. That implies that all sums over intermediate states run over some 
finite set of one-electron states. The latter should effectively account for both discrete 
and continuous part of the spectrum of Dirac-Fock operator (||) . 

In this calculation basis set included 167 four-component orbitals: 1 — 23s, 2 — 22p, 
3 — 22<i, 4 — 15/, 5 — 15g and 6 — 13h. Core orbitals as well as valence orbitals 
6s, 7s, 6p, 7p, and 5d, 6d were obtained by solving Dirac-Fock equations for appropriate 
configurations, while higher virtual orbitals were formed, as described in |2l], [22), [I], |J. 
After that operator (||) was diagonalized in order to have quasi-Dirack-Fock basis set. 

This basis set was used to solve RPA equations for HFS operators and to calculate 
diagrams for effective Hamiltonian and effective HFS operators. After that, the 
solution of Equation (Q) was found in the CI approximation. On this stage it was 
possible to truncate our basis set to 7 — 15s, 6 — 15p, 5 — 14<i, 5 — 10/ and 5 — 7g 
orbitals. We made complete CI on this shorter basis set. 

We calculated effective Hamiltonian _ff c ff in the second order of MBPT. That 
means, that zero order approximation ( p"4| ) for the Green's function was used in 
(||). The choice of the constant W suggested by Equation (|l0|) resulted in some 
underestimation of the binding energy for levels of the configuration 5d6s and for l P° 
level of the configuration 6s6p. For these levels W = Wo + 0.4 au was used. Our 
final results for the ground state are given in table |l| and for several low-lying levels 
in table g 

On the next stage we used corresponding wave functions to calculate HFS 
constants. These calculations are similar to those of the papers Q and || . Results are 
given in tables || and [|. In analysis of the theoretical accuracy, which is done in the next 
section, it is very important to know the scale of different contributions to the final 
answer. Therefore, in tables || and || we give results of the Dirac-Fock calculation, of 
the conventional two-electron CI and of the two-electron CI with effective Hamiltonian. 
Final values include effective operator corrections ( |l7j ) and normalization correction 
(|U). The latter usually decreases the answer by l%-2%. 

4. Discussion 

We have pointed out above, that the effective Hamiltonian can be safely used only 
for the energy levels below the core excitation threshold. For Yb this threshold lies 
at 23189 cm^ 1 above the ground state and many of the levels from the table |^ lie 
higher. Nevertheless, theoretical spectrum is quite good up to the level 1 PJ(6s7p), 
which appears 1634 cm -1 below its experimental position. This huge discrepancy can 
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be attributed to interaction with the J = 1 level at 38422 cm -1 , which is supposed to 
belong to the configuration 4/ 13 5d 2 6s Q. This level can also interact with 3 PJ(6s7p), 
what makes calculations for this level less reliable. For other levels from the table || 
there are no close levels with the same J and the same parity, which correspond to 
the excitations from the 4/-shell. 

The perturber levels, discussed above, correspond to the poles of the Green's 
function Rq{E) and, thus, to the poles of the operator S(-E). As long as we use 
approximate Green's function ( p"4| ) instead of the exact one, the poles of X(-E) are 
shifted from the physical poles. In order to use effective Hamiltonian, we have to 
keep far enough from the poles of the exact Green's function and from the poles of 
the approximate Green's function. As we go further up above the core excitation 
threshold, both sets of poles become denser. For this reason, it is hardly possible to 
use effective operator technique there. 

Let us now proceed to the calculations of HFS constants presented in the tables || 
and |J. In these tables the first two rows correspond to the non-perturbative part 
of the calculation, while two lower rows include MPBT corrections associated with 
the effective Hamiltonian and the effective HFS operators. One can see, that average 
MBPT contribution is about 30% for magnetic constant A and about 40% for electric 
quadrupole constant B. Thus, theoretical accuracy for the constant A is normally 
higher. Indeed, we have neglected higher order corrections to the effective Hamiltonian 
and to the effective HFS operators, that can be justified only if MBPT corrections are 
small. 

For 6 levels ( 3 D 2 (5<26s), 3 D 3 (5<26s), 3 Si (6s7s), 3 Di (6s6d), 3 D 3 (6s6<2), and 
3 P2 (6s7p)) MBPT corrections to the constant A are less than 25%. For all these 
levels the difference between the theory and experiment is within 3%. For the constant 
B the smallest MBPT corrections (about 35%) correspond to levels 3 Di (6s6d) and 
3 P° (6s6p), where we have agreement with the experiment within 2%. Constant B for 
3 Si (6s7s)- state differs from zero only because of the configuration interaction. 

On the other hand, when MBPT corrections are 40% or more, the accuracy of 
the theory becomes uncertain. The most striking discrepancy with the experiment 
takes place for levels X D 2 (5d6s), 3 D 2 (6s6d) and : Pf (6s6p). For all of them MBPT 
corrections appear to be 40%-50%. 

The only exception of this rule is the level 3 PJ (6s7p), where MBPT corrections 
for both constants are about 30%, but even the sign of the constant B differs from 
the experiment. As we pointed out above, this level can interact with the perturber 
level 4/ 13 5d 2 6s(J = 1), which is located at 38422 cm -1 . Constant B for this level is 
relatively small and even small admixture of the perturber can result in the change of 
its sign, while the larger constant A should be less affected by such a mixing. 

5. Conclusion 

In this paper we have checked applicability of the CI+MBPT method for ytterbium. 
The first core excitation here lies only 23189 cm -1 above the ground state, which 
seems to restrict us to only few lowest levels. Actually, the method works quite well 
up to 38930 cm -1 , where the first significantly perturbed level arise. 

When the method is used for calculations of HFS constants, the accuracy of 
calculations depends on the scale of MBPT corrections. When total MBPT correction 
contributes less than 25% to the value of the constant, the accuracy appears to be 
better than 3%. On the contrary, when MBPT corrections account for more than 
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40% of the theoretical value, the accuracy becomes uncertain. In this case one can 
use MBPT correction only as a rough estimate of the theoretical error. We expect, 
that similar relation holds for other one-electron operators as well. 

For the studied levels of Yb the average MBPT correction is approximately 30% 
for the constant A and 40% for the constant B. Therefore, precision calculations of 
the constant A are possible for the majority of levels, while precision calculations of 
the constant B are possible only as exception. 
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Table 1. Two-electron energies E v& \ of the Yb ground state in different 
approximations (au). MBPT corrections were calculated with W = Wo (see 
Equation ((Hp). 



DF a 


MBPT b 


CI C 


CI+MBPT d 


Experiment 


0.606810 


0.654308 


0.632398 


0.677601 


0.677584 



a Single-configuration approximation. 

b Single-configuration approximation with MBPT correlation. 
Conventional CI method for two outer electrons. 

d CI+MBPT method: two-electron CI with effective Hamiltonian formed within 
second order of MBPT. 

e This is the sum of the first two ionization potentials of Yb |23[. 



Table 2. The low- lying levels of Yb in V N approximation. -E va i is the two- 
electron energy and A is the interval from the ground state. The multiplet 
splittings are given in parentheses. Levels marked with f were calculated with 
W = Wo + 0.4 au and the levels marked with J were calculated with W = 
Wo + 0.1 au. 

Theory Exper. a 



Level 


E val (au) 


A (cm -1 ) 


A (cm" 1 ) 


Even states 


3 Di (5(26s)t 


-0.566086 


24441 


24489 


3 D 2 (5<f6s)f 


-0.564919 


24697 (256) 


24752 (263) 


3 D 3 (5<26s)t 


-0.562415 


25247 (550) 


25271 (519) 


X D 2 (5(26s)t 


-0.551500 


27642 


27678 


3 Si (6s7s) 


-0.529514 


32501 


32695 


1 S (6s7s) 


-0.521725 


34210 


34351 


3 Di (6d6s) 


-0.496053 


39845 


39809 


3 D 2 (6d6s) 


-0.495902 


39877 (33) 


39838 (29) 


3 D 3 (6d6s) 


-0.495384 


39991 (114) 


39966 (128) 


X D 2 (6d6s) 


-0.494733 


40135 


40062 



Odd states 



3po 


(6s6p)t 


-0.598858 


17282 


17288 


r i 


(6s6p)| 


-0.595599 


17997 (715) 


17992 (704) 


3po 
i 2 


(6s6p)t 


-0.587613 


19750 (1753) 


19710 (1718) 


Ipo 


(6s6p)f 


-0.563354 


25074 


25068 




(6s7p) 


-0.504463 


38000 


38091 


3pg 
, 1 


(6s7p) 


-0.504096 


38080 (80) 


38174 (83) 


3po 

i 2 


(6s7p) 


-0.502388 


38455 (375) 


38552 (378) 


Ipo 

r 1 


(6s7p) 


-0.500224 


38930 


40564 



a Rcf. H 
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Table 3. Magnetic dipole (A) and electric quadrupole (B) hyperfine structure 
constants of low-lying even-parity levels for 173 Yb. The electric quadrupole 
moment is taken to be 2.80 b. 





3 D 1 


3 D 2 


3D 3 


X D 2 


3 S 1 


3 D 1 


3 D 2 


3 D 3 


X D 2 




(5d6s) 


(5d6s) 


(5d6s) 


(5d6s) 


(6s7s) 


(6s6d) 


(6s6d) 


(6s6d) 


(6s6d) 


A (MHz) 


DF 


447 


-200 


-348 


-18 


-1225 


490 


-206 


-337 


28 


CI 


443 


-288 


-348 


62 


-1489 


633 


-468 


-437 


238 


Heft 


588 


-398 


-469 


105 


-1910 


830 


-677 


-567 


380 


Total 


596 


-351 


-420 


131 


-1938 


838 


-683 


-569 


392 


Exper. a 


563(1) 


-362(2) 


-430(1) 


100(18) 


-1879(10) 


818.7(4) 


-732.5(4) 


-559.9(5) 


438.5(4) 



B (MHz) 



DF 


151 


206 


353 


379 





32 


15 


76 


81 


CI 


156 


229 


368 


676 


0.3 


38 


55 


92 


115 




219 


314 


492 


693 


0.2 


38 


62 


92 


115 


Total 


290 


440 


728 


1086 


0.2 


58 


96 


150 


184 


Exper. a 


335(1) 


482(22) 


909(29) 


1115(89) 


<3(18) 


59.3(2) 


52.5(2) 


139.6(3) 


142.2(2) 



a HFS constants for configuration 5<i6s were measured in [g4J, for 3 Si-state in |g5J, and for configu- 
ration 6s6d in [f!6| 



Table 4. Magnetic dipole (A) and electric quadrupole (B) hyperfine structure 
constants of low-lying odd-parity levels for 173 Yb. 





3 po 


3 po 


1 po 


3 po 


3 po 




(6s6p) 


(6s6p) 


(6s6p) 


(6s7p) 


(6s7p) 


A (MHz) 


DF 


-664 


-527 


- 25 


-667 


-507 


CI 


-765 


-556 


98 


-1044 


-666 


H eB 


-1075 


-741 


187 


-1468 


-862 


Total 


-1094 


-745 


191 


-1488 


-871 


Exper. 


-1094.0(7) a 


-738 b 


59 b 


-1144(l) c 


-854(2) c 




-1094.2(6) d 










B (MHz) 


DF 


-454 


715 


811 


-90 


150 


CI 


-533 


860 


428 


-80 


163 


H eB 


-633 


1013 


642 


-82 


192 


Total 


-822 


1335 


848 


-111 


264 


Exper. 


- 826.5(l) a 


1312 b 


605 b 


12(4) c 


267(25) c 




-827.2(5) d 











a Rcf. [Q; b Rcf. H; c Ref. d Ref. || 



